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Abstract: 

Hepatitis C virus (HCV) is a source of chronic hepatitis C infection that can destroy the liver 
when proper therapeutic intervention is not executed early. For the purpose of studying 
the dynamics of the disease during acute infection, a deterministic mathematical model 
incorporating the effect of immune system and therapy has been formulated. The model is 
a system of non-linear ordinary differential equations. The disease-free equilibrium point 
E, during therapy has been computed. Furthermore, the basic reproductive number R,and 


effective reproductive number R, have been computed using the next generation operator 
method. The sensitivity indices of R, relating to each parameter in the model were 
computed. It is found that R,is most sensitive to the drug inhibiting efficacy of viral 


replication parameter € . Thus, it implies that an increase in the value of £ will reduce the 
rate of viral replication leading to the decrease or eradication of the HCV morbidity. By 
using the estimated parametric values, we found that R,<R,, which attests that 


therapeutic strategy is absolutely effective in reducing intensity of the disease and hence 
preventing evolution to chronic state. 


1. Introduction 

HCV is a small positive-stranded RNA microbe belonging to the Hepacivirus genus in the 
family Flaviviridae (Taxonomy, 2005), which mutates so rapidly that there is no vaccine 
(Colina et a/, 1999). It is the source of hepatitis C virus infection (HCVI) that can destroy 
the liver if no proper therapeutic strategy is implemented. It is approximated that about 
130 to 170 million people are infected with HCV worldwide (Lavanchy, 2009). Usually, 
HCVI is categorized as acute hepatitis C infection (AHCI) and chronic hepatitis C infection 
(CHCI) depending on the disease evolution. The AHCI signifies the first 6 months of 
infection following acquisition of HCV whereas CHCI refers to a period subsequent to the 
AHCI phase, which is marked by the persistence of HCV RNA in the blood. This is 
irrespective of whether patients are symptomatic or asymptomatic during AHCI. It is also 
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evident that around 20% to 50% patients with AHCI spontaneously clear the virus within 6 
months (Kamal, 2008) while others progress to CHCI. It implies that more than 50% of 
AHCI cases undergo CHCI, which can be a source of evolution of cirrhosis and 
hepatocellular carcinoma to some of them. Many people with AHCI are asymptomatic, 
which makes the diagnosis of the disease become difficult while symptomatic patients most 
commonly show fatigue. 


The channels of HCV transmission are primarily dependent on exposure to infected blood, 
which incorporate vertical transmission (Owusu-Ofori et a/, 2005), needle stick injuries 
(Xia et al, 2008), intravenous drug use (IDU) (Tohme and Holmberg, 2010), sexual 
behaviors(VJafari etal, 2010) and body piercings (Lam et a/, 2010). In developing countries 
(DCs), HCV is transmitted mainly through contact with infected blood and blood products 
in healthcare provisions centers. In most developed countries, it results from IDU; it can be 
a risk source of the transmission in the DCs as well (Aceijas and Rhodes, 2007;Nelson et a/, 
2011). 


We know that numerous antivirals have been clinically tested for therapy of HCV. The trials 
have revealed that therapy of CHCI with different antivirals results into different SVR rates. 
For instance, it was formerly treated using conventional interferon alpha that led to the 
SVR rate of at least 20% of patients after 48 weeks of therapy, which later was improved to 
the SVR of 38% to 43% of patients using interferon alpha-2b plus ribavirin combination 
therapy (McHutchison et al, 1998; Poynard et al, 1998). Combination therapy with 
peginterferon alpha 2a or 2b plus ribavirin combination therapy further improved the SVR 
rate to approximately 50% of cases (Fried et a/, 2002; Muir et al, 2004; Kanda etal, 2010), 
which has been the standard therapy for CHCI. Unlike AHCI, chronic hepatitis c is 
associated with a worsening prognosis, which requires more demanding therapy and 
longer therapy period leading to intolerability owing to developing side effects. Conversely, 
drugs prescribed for therapy of CHCI can also be assigned to people with AHCVI although 
interferon-based monotherapy is just sufficient to produce SVR of more than 50% of cases 
(Jaeckel et al, 2001; Kamal et al, 2006). Monotherapy is far better tolerated due to fewer 
side effects, relatively cheaper, more suitable and takes shorter time than combination 
therapies for chronic CHCI (Myers etal, 2001).Hence, therapy of patients with HCV during 
acute period of infection significantly reduces side effects and disease development to 
chronic hepatitis(Alberti et a/, 2002). 


Although this paper presents a deterministic mathematical model for the effect of therapy 
on dynamics of HCV in vivo during acute phase of infection, patients with AHCI are diverse. 
It is very difficult to identify patients who require therapy owing to varying rates of 
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spontaneous clearance of the viruses and the presence of asymptomatic cases (Busch and 
Shafer, 2005). Spontaneous clearance of the virus can take place either during AHCI (Kamal, 
2008) or beyond this period (Grebely et al, 2014).Also, some patients are entirely 
asymptomatic having normal or slightly elevated serum alanine aminotransferase (ALT) 
levels despite HCV-RNA positivity after known exposure. Others show acute illness and 
have high serum ALT levels and noticeable jaundice. Thus, we have considered these 
challenges in the development of the model, i.e. a patient can either clear the virus 
spontaneously or not during acute infection. 


Certainly, mathematical modeling of viral dynamics with therapy provides us with 
understanding and quantifying biological processes that govern changes in the viral load 
and accompanying biomarkers such as ALT levels before and after treatment (Ribeiro et al, 
2003).For our case, mathematical models can help to provide answers to biological 
questions regarding pathogenesis, the dynamics of HCV, immune response and efficacy of 
drug therapy. Furthermore, we see that modeling HCV dynamics with therapy has 
deepened our understanding of the virus pathogenesis and guided drug development 
(Chatterjee et a/, 2012). Various models have been formulated by several researchers in 
order to understand the efficacy of drug therapy on the dynamics of HCV. Thus, Neumann 
et al(1998) formulated a simple model to investigate the dynamics of HCV during therapy 
by adopting a model of human immunodeficiency viral infection (Wei et a/, 1995; Perelson 
et al, 1996). The model had a system of three nonlinear differential equations which 
represented a constant population of target hepatocytes, infected hepatocytes and the virus, 
and it was used to approximate the rates of viral clearance and loss of infected hepatocytes 
by fitting to the model the decline of HCV RNA level within infected people in the first 14 
days of therapy. Since then modeling of the virus dynamics has been very important in 
examining the HCV RNA level decay during therapy, which is very well expounded by 
Perelson (2002) and Perelson et a/.(2005). However, it could not include some observed 
HCV RNA dynamic profiles under therapy, i.e., it assumed that there is no proliferation of 
target and infected cells and no spontaneous recovery of infected hepatocytes. The 
extended model of Dahari et a/(2007) incorporated proliferation of target and infected 
hepatocytes without spontaneous recovery of the infected hepatocytes. A model that 
includes proliferation and the spontaneous recovery of infected hepatocytes was used to 
investigate the dynamics of HCV for primary infection in chimpanzees (Dahari et al, 
2005).Nevertheless, these models do not include the role of immune response on the 
dynamics of HCV. 


The Neumann et a/(1998) model and extended model of Dahari et a/(2007) assumed that 
the drugs had two major effects on the dynamics of HCV: reduction of infections and 
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blockage of virus replication. We observe that the drug blocking efficacy of virus replication 
increases with drug dosage (Neumann et a/, 1998), where at large values of blocking 
efficacies the drug does not directly reduce infections i.e. the drug efficacy in reducing 
infections becomes negligible. Thus, it is required to consider only blockage of the virus 
replication. In this paper, we extend the model of Ismail eta/ (2016) to investigate the 
effect of interferon alpha-2b monotherapy on the dynamics of HCV by adopting therapy 
regimen of Neumann et a/(1998). 


2. Model Formulation 

In this section, we formulate a deterministic model incorporating three different 
populations: hepatocytes population, Hepatitis C virus population and the CD8+ T cells 
population. But, we classify the hepatocytes population as susceptible and infected classes 
whereby the susceptible class is a group of hepatocytes (liver cells) exposed to the HCV 
infection, and are equally likely to be infected, while the infected class is a group of 
hepatocytes infected with HCV. Moreover, the HCV population refers to a group of viruses 
that is responsible for the infection whereas the CD8+ T cells population is a collection of T- 
Lymphocytes of the immune system that is liable for the destruction of the infected 
hepatocytes. We have extended the model proposed by Ismail, et a/ (2016) to include the 
effect of antiviral monotherapy. 


In order to achieve our goal, we make the following assumptions: 

(1) The model is formulated to describe the dynamics of HCV and immune system during 
acute phase of infection; (2) susceptible hepatocytes are equally likely to be infected by the 
viruses and infected hepatocytes; (3) the susceptible hepatocytes are generated by 
proliferation of existing hepatocytes and immigration (differentiation of liver cells 
progenitor or bone marrow cell) at a constant rate; (4) these susceptible hepatocytes and 
infected hepatocytes have the same natural mortality rates; (5) the virions and the CD8+ T 
cells have different natural death rates; (6) the CD8+ T cells destroy the infected 
hepatocytes at a constant rate and die at a constant rate due to infection; (7) the CD8+ T 
cells are produced at a constant rate and (8)we assume that the person with HCV infection 
receives highly dosed interferon alpha-2b monotherapy that blocks viral replication within 
infected hepatocytes. (9) Also, we assume that the person with HCV infection can either 
clear the virus spontaneously by a noncytolytic process or not during therapy. 

All state variables and parameters are positive; they are introduced and concisely 
described in Table 1 and Table 2 respectively. 
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Table 1: State variables and their descriptions 


Variable Description 
S(t) The susceptible hepatocytes at time t 
I(t) The infected hepatocytes at time t 
V@) The hepatitis c viruses at time t 
T(t) The CD8* T cells at time t 


For brevity, we henceforth denote the susceptible hepatocytes sub-population by S; the 
infected sub-population by Z ; HCV viral population by V and the CD8+ T cells population by 
T .The state variables and their descriptions are summarized in Table 1. 


Table 2: Parameters and their descriptions 
Parameter Description 
Qo Per capita infection rate 
Per capita production rate of viruses from the infected hepatocytes 
Rate at which the CD8+ T cells destroy the infected hepatocytes 
Per capita production rate of susceptible hepatocytes 
Per capita production rate of the CDB+ T cells 
Per capita natural death rate of susceptible and infected hepatocytes 
Per capita natural death rate of viruses 
Per capita natural death rate of CD8+ T cells 
Rate of spontaneous cure of infected hepatocytes by a noncytolytic 
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process 
Per capita death rate infected hepatocytes due to HCV infection 
E Fraction by which antiviral drug reduces viral production rate 
Eez Maximum CD8+ T cells population level 

Healthy hepatocytes Sare constantly produced at the rate II and die naturally at a 
constant rate d; the S cells are infected at the rate that is proportional to the product 
SV ,with a constant of proportionality œ.The infected hepatocytes 7 have a constant 
natural mortality rate of d .The hepatitis C viruses V are produced from these infected cells 
at a constant rate o viruses per infected cell per day, once produced they die naturally at 
the rate of c.Then the CD8* T cells destroy the infected cells 7 at the rate proportional to 
the product 7 and T, with a constant of proportionality 8. In the presence of HCV, the 
CD8+ T cells are instigated by the CD4+ T cells and constantly supplied at the rate g,and 


die naturally at a constant rate b . When the person infected with HCV receives interferon 
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alpha-2b monotherapy, the drug inhibits virus production from infected hepatocytes 7 bya 
fraction €. 


2.1 Model Compartment and Dynamics 

If we consider the assumptions and the symbols representing the parameters and state 
variables from Table 1 and Table 2 respectively, we can illustrate the description of the 
model dynamics in the form of a compartmental diagram as in Figure 1. 





gV(-T/T_) (te)ol 


Figure 1: Compartmental diagram for the hepatitis C virus model with immune system 
including drug therapy. 
2.2 Equations of the model 
Based on the assumptions made and relationships that exist between the state variables 
shown in Figure 1, we formulate a system of four non-linear ordinary differential equations 
describing the dynamics of hepatitis c virus in the presence of drug therapy. 


dS 


a a E (1.1) 
dI 

pore ee ee (1.2) 
dV 

—=(l-e)ol-cV 1.3 
F (l-e)ol —c (1.3) 
dT T 

— = eV(1-——) -bT 1.4 
7 gV( 7 ) (1.4) 


max 


where the initial conditions are S(0)> 0,7(0)> 0, v(0)> 0and 7(0)>0 
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2.3 Basic properties of the model 
2.3.1 Invariant regions 
Since the system of equations (1.1)-(1.4) describes modeling of susceptible hepatocytes, 
infected hepatocytes, HCV and the CD8+ T cells with the effect of drug therapy, we assumed 
that all state variables and parameters used in the model are non-negative vt > 0. 
According to Figure 1, we observe that the hepatocytes population has two classes S(t)and 
I(t) which can be combined together to form N(t), which denotes the total hepatocytes 
population. That is, 
N(t)=S(t)+1@) (2) 

Then we can determine the appropriate feasible region where all state variables are non- 
negative through the following theorem: 
Theorem 1: A// forward solutions of the system (1) are contained in the region® € R*, 
Vt2Oand®=T, OT, Ul, € R? xR! x Ri, where 

T, =(S,1)eR?:S4+1<N 


T, =Ver! 
I, =TeR! 
and © is the positive invariant region for the whole system. 


Proof: 
We prove the theorem by initially determining the invariant region within which the 
solution for each population are feasible Vr > 0. 


Hepatocytes population 
We have to determine the invariant region I’, of the system (Hepatocytes) containing the 
feasible solutions Vr > 0. 


Let T, = (S, I) e R? be any solution of the system with non-negative initial conditions. 


From (2), we have: 


aN = ds + dI (3) 
dt dt dt 
Substituting (1.1) and (1.2) into (3) produces 
AA I+ qI- oSV - dS] + [08V ~ AIT -dl — pd - q1) 
=M- Nd - PIT — Wl (4) 
From (4), we have: 
dN 


—<II-Nd 
dt 
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DN Net (5) 
dt 


The general solution of the first order ordinary differential inequality (5) is: 


gi 6 
NOS or ATC (6) 


where N, is the initial hepatocytes population size. 


From (6), we deduce that 
II 
N(t)s marto (7) 


Thus, the feasible solutions for the hepatocytes population in the system (1) are positively 
invariant in the region: 


T, =(N@:N(@)<max{Ny,—)} 


HCV population 
We have to determine the invariant region I), of the system (HCV) containing feasible 


solutions Vt >0. Let T, =V e R} be any solution with non-negative initial conditions. 


Since we know that 
From (2) and (7), we deduce that 
II II 
I(t) < max(N,, ru >I(t)< Fi (8) 
Substitution of (8) into (1.3) produces 
dV Z (-e)oil _ 
dt d 


Lg teal 
dt 


cV 





(9) 


The general solution of the first order ordinary differential inequality (9) is 


(l—«)oi1 (l-«) 
cd 


V(t) < +(V, — Jet" exper) (10) 
c 


where V, is the initial viral population size. 

From (10), we deduce that 

d-e)d I 
cd 


Thus, the feasible solutions for the viral population in the system (1) are positively 


V(t) < max{V,, ) (11) 


invariant in the region: 
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(-es)dl 


Ty, ={V(t): V(t) < max{V,, 
cd 


k 

CD8+ T cells population 

We have to determine the invariant region I’, of the system (CD8+ T cells) containing 
feasible solutions vt > 0. 


Let T, =T e R! be any solution with non-negative initial condition 


From (11), we deduce that 














V(t) < Q-ə)dl (12) 
cd 
Substitution of (12) into (1.4) produces: 
dT < g(d-e)dll a T )-bT 
dt cd T 
T AL L (13) 
dt cdT ax cd 
The general solution of the first order ordinary differential inequality (13) is 
MOE EASRA Ty A E E Hy) (14) 
g(—-e)ol1+bcdT,,,, g(l—e)ol1+bcdT,,,, cdT ax 
From (14), we deduce that 
1-e)oll 
T(t) < max{T), sts) } 





g(l—e)ol1+bcdT,,,. 
Thus, the feasible solutions for the CD8+ T population in the system (1) are positively 
invariant in the region: 
g(l-e)ol1 


ne AT SO Fae ean De 





}} 
Thus © =T, Ul, UT, € R? x R! x Ri, for 
II 
p UNAS RARR, 14 >0}, 


d-e)oiI 

ae }} 
gd-e)dl 

g(l—e)ol1+bcdT.,,. 


From this, we conclude that the model system (1) is positively invariant in the region 


Iy ={V(t): V(t) < max{V,, 


I; ={T@:T@ <max{Th, 





i 


® .Hence it is epidemiologically and mathematically realistic. 


2.3.2 Positivity of solutions 
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Since the system (1) refers to modeling of populations, where all state variables and 
parameters are assumed to be non-negative Vt > 0, we have to solve each equation in the 
system to test for positivity. We achieve this thorough the following theorem: 
Theorem 2: /f the initial values ofa given system are {S(0),1(0),V(0),T(0) € R *}20 then 
the solution set {(S(t),1(t),V(t),T@)} consists of positive entitiesVt 20. 
Proof: 
We test for positivity of each state variable. 
From (1.3), we have: 

ay 2—cV 

dt 

dV 

vv = —cdt (15) 
Integration of (15) produces 

ln V > —ct+ K 
That is, V(t) 2V,, exp(—ct) (16) 


where V, is the initial viral population size. 


From (16), we have: 
Att =0,exp[—c(0)] =1.S0,V(O) = V, 2 0.1f t >0,exp[—c(4)] > O since c > 0.Thus, we find that 


V(t) = V exp(—ct) = 0, Vr >0 


From (1.1), we have: 


1 >-(@V +d)S 
dt 


o> (oV +d)dt (17) 


Integration of (17) produces: 
InS >-(@V + d)t+ K 


That is, S(t) 2 S, exp[ (V(t) + d)t] (18) 
From (18), we have: 

Att =0,exp[-(@V (0) + d)(0)]=1.So0, S(O) = S, 2 0.If t > 0, then exp[—(@V (t) + d)t] > O since 
ov (t)+d >0. Thus, we find that S(t) > S, expl (oV + d)t] 20, Vr>0 


From (1.4), we have: 


dT, (EV HOT 
dt T ax 
dT , (BV + Tw 
T T, 


max 





~\T 





‘dt (19) 
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Integration of (19) produces: 
gV +bT ax 


InT >-( \t+K 
That is, 

T(t) = Ty exp ET (20) 
From (20), we have: 
Att =0, E 1 Fax y(Q =1.So, T(0) =T, = 0.1ft >O then exp oY Plies s] > Osince 


gV+DT,, gV+DT., 


< > 0.Thus, we find that T(t) 2 T, exp[—( ~)]=0,Vvr=0 
From (1.2), we have: 
a >-(PT+d+q)I 
dt 
dI 
A EEU EUA (21) 
Integration of (21) produces: 
InJ>-(PT+d+q)t+K 
That is, 
I(t)> 1, exp HET +d +q)t] (22) 


From (22), we have: 

Att =0,exp[~ 8T (0) + d + q)(0)] =1.S0, I(0) = J, 2 0.Ift >0, then exp[—(G7 (1) + d +q)t]>0 
since £T(t)+d+q>0 Thus, we find that Z(t) = I expl (£T +d +q)]20, Vr=0, 

Since the solution set {S(z),/(t),V(t),7(t)} consists of positive entities Vr >0, we conclude 


that the model system (1) is epidemiologically and mathematically realistic(Hethcote, 
2000). 


3. Model Analysis 

In this section, we determine the disease free equilibrium point of the system (1) and the 
basic reproductive number in the absence of therapy. We also determine the effective 
reproductive number for the purpose of examining the impact of antiviral therapy on the 
HCV morbidity. 


3.1 Disease Free Equilibrium Point (DFE) 
We determine the disease free equilibrium point by setting the equations of the system (1) 
equal to zero. 
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I1+ ql —a@SV —dS =0 (23.1) 
OSV — ZIT — dI — pl — ql =0 (23.2) 
(l—e)ol -cV =0 (23.3) 
sV(l-)-bT =0 (23.4) 
From equations (23.1), (23.3) and (23.4), we obtain 
. H+qř i 3 T aV” 
oe Le iS TENE pita. (24) 
oV“ +d (l-e)o gV’ +bT ax 
At the disease free equilibrium point, V* =0.Then, from (24), we obtain 
T x (0 
œo(0)+d d (l-é)o g0) +bT ax 
Hence the disease free equilibrium point of the system (1) exists and is given by 
II 
Es = F ,0,0,0) (26) 


According to (26), the absence of hepatitis c virus implies absence of infected hepatocytes 
and therefore the CD8+ T cells are not activated to destroy the infected hepatocytes. 


3.2 The Basic Reproduction Number, R, 
Definition 1:7he basic reproduction number ,R,,is defined as the average number of 


secondary infections produced when an infected individual is introduced into a host 
population where all individuals are susceptible in the period of infection (Dietz, 1975; 
Diekmann et al.; 1990, Van den Driessche and Watmough, 2002). 


The basic reproduction number is a very important parameter as it can be used to 
determine whether the disease prevails in the community or dies out. If R, <1, the average 


number of new infected individuals produced by an infectious individual is less than one 
during its infectious period, implying that the infection cannot grow. Conversely, R, >1 


implies that the average number of new infections produced by an infectious individual 
greater than one, and therefore the disease spreads in the population. It is also very 
important in the process of analyzing sensitive parameters which drive the dynamics of the 
disease and stability analysis of the disease free equilibrium. 


To determine the basic reproduction number, we find it important to recognize new 
infections from all other changes in the population. To achieve this, we apply the next 
generation operator proposed by Van den Driessche and Watmough (2002). Then we 
consider the model system (1) without therapy, i.e.g =0. 
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Then in the absence of therapy, we have the system (1) reduced to the system (27) 





E tae as 

dt 

E pid 21 

dt 

dV ae. 
—=ol-—cV 

dt 

dT T 

“= gV(l-——_) -bT 

zi gV( T ) 


max 
Let f,(x) be the rate of appearance of new infection in compartmenti,v, (x) be the rate of 
transfer of individuals out of compartmenti and v,' (x) be the rate of transfer of individuals 
into of compartmenti by all other means. Furthermore, we assume that each function is 
continuously differentiable at least twice in the region ®.The HCV model system comprises 
nonnegative initial conditions with the following equations: 

x, =F (x = fi- y; (x), i =1,2,...,7n 


where v, =v; -v` 
Then we consider equations (1.2) and (1.3) of the system (1) to determine the expressions 
for f, and v,as follows: 
oSV (PT +d+yu+q)l 
le] | 
Let F be a non-negative matrixnxn and Y bea non-singular N -matrix 
At the disease free equilibrium point E,, we find that 





ex OREN shor 
Ox c o0 
Ox, 0 c 
Then, we have 
ot 
y= d+u+q 
0 1 
c 


We determine the next generation matrix FV ' developed by Castillo-Chavez et a/.(2002) 


as 
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0 all 
FY" = cd 
o 
— 0 
d+u+q 


By Van den Driessche and Watmough (2002), the basic reproductive number is a spectral 
radius of the next generation matrix FY~'.Henceforth, to determine the basic reproductive 
number R,, we compute the spectral radius (p) of this next generation matrix from 

R, = P(FY*)=max(/,,/,), 
where 4 and /,are the eigenvalues of the next generation matrix FY~™' 


Thus, the eigenvalues, lambda, of Fy“ are obtained by 


[FY —Al|=0 
_, al 
cd =0 
That is, ro -A 
d+u+q 
> aoll 





=0 
cd(d + u+q) 


Hence the basic reproductive number is given by 


R,= (ee 0 (28) 
cd(d+ u+q) 


3.3 The Effective Reproduction Number, R, 
Definition 2:7he effective reproduction number, R,,is the average number of infections 


caused by a single infectious individual introduced in a community where intervention 
strategies (in our case therapy) are administered (Okuonghae and Korobeinikov, 2007; 
Okuonghae and Aihie, 2008). 

We compute the effective reproduction number R,of the system (1) by applying the 


method used for R,,where R, is the spectral radius of FY ` denoted by R, = p(FY ').Thus 


the effective reproduction number is given by 


jee E a (29) 
cd(d+u+q) 


R, =/0-8)R, (30) 


From (29), we have: 
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The drug blocking efficacy, € , of viral replication can take any value in the interval O<¢<1, 
whereby its role can be discussed in three sub-intervals: ¢=0,0<é<landé=1.When 
£ =0, the result in (30) reduces to the form R, = R,,which implies that drug therapy has 


no impact on the HCV morbidity. Conversely, ife =1, it becomes R, =0, implying that the 


disease vanishes due to therapy, i.e. the drug is 100% effective in exterminating the disease. 
But, any value of <€ taken from the sub-interval 0<é<1 changes it into R, < R,.This 


reflects therapy efficacy to reduce the impact of the disease. 


4. Numerical Sensitivity Analysis 

Mathematical epidemiology is a study that helps in understanding the dynamics of an 
infection in order to control it by focusing on some sensitive parameters. This can be 
achieved by performing sensitivity analysis based on the model parameters. The sensitivity 
analysis describes how each parameter influences R, .It helps to identify the most sensitive 


parameters. With small variations in the values of the parameters, we can identify the 
parameters that highly influence R, whereby proper control measures can be taken(Chitnis 


et al, 2008). We accomplish this by computing the sensitivity indices of R, relating to all 


parameters in it. In this case, the estimated values of the parameters employed were 
adopted from some literatures and others were just estimated, which are altogether 
itemized in Table 3. 


Table 3: Estimated values of parameters used to compute sensitivity indices of R, 


Parameter Value Unit Source 
b 0.02 day! Avendano et a/(2002) 
c 10 day! Estimated 
d 0.00014 day! Estimated 
& 0.0003 day! Avendano et a/(2002) 
q 2 day! Estimated 
2 0.00000001 virus cell-tday-t Dahari et al. (2005) 
u 0.486 day Estimated 
E 0.96 Neumann et a/(1998) 
o 6 virus cell-tday-t Estimated 
o 0.00001 virustmltdayt Estimated 
TI 100 cells ml-!day! Estimated 
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The normalized forward sensitivity index refers to the ratio of relative change of a variable 
to the relative change in a parameter. If the variable is a differentiable function of the 
parameter then the sensitivity index is defined as follows: 


Definition 3: The normalized forward sensitivity index of variable p that depends on 


parameter q is defined as: 





õp 4 
P — x 31 
Xa a p (31) 
Replacing p in (31) by R,, we obtain 
OR, q 
Re = e p aea 32 
Aa ôq R, ea 


According to (32), the sensitivity indices of R, are computed with regard to all parameters 
in it. This is achieved by replacing q by a parameter to obtain 
OR OR 
ye =—x kae —12.0000, y} = — x ce —0.5000 and so on 
Oe R Od R 


The sensitivity indices of R, pertaining to all parameters are computed using the same 





approach and then summarized in Table 4. 


Table 4: Sensitivity indices of R, relating to all parameters 


Parameter Sensitivity index 
E —12.0000 
d — 0.5000 
c —0.5000 
o + 0.5000 
oO + 0.5000 
II + 0.5000 
q — 0.4022 
u — 0.0977 


According to the sensitivity indices displayed in Table 4, we observe that the parameters 
for infection rateq@, the viral replication rate o and rate of recruitment of susceptible 
hepatocytes II are the most positively sensitive parameters. This implies that increasing 
the values of these parameters will increases R, and vice versa. For example, increasing the 


value of aby 20% will increase R, by 10% while decreasing w@by 20% will decrease R, by 
10%. 
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The drug blocking efficacy of virus replication € , the natural mortality rate of susceptible 
and infected hepatocytes d and the natural death rate of the virusc are the most negatively 
sensitive parameters. This implies that increasing the values of these parameters will 
decrease R, and vice versa. For example, increasing ¢ by 10% will decrease R, by 120%. 


Also, we observe that the rate of spontaneous recovery of infected hepatocytes g and the 


death rate of infected hepatocytes due to infection w are less sensitive negative parameters. 


5. Numerical Simulations and Discussion 

In this section we present the numerical simulations of the system (1).We performed 
simulations based on the estimated parametric values itemized in Table 3 and this was 
accomplished through the application of ode45 MATLAB's standard solver for ordinary 
differential equations (ODEs), which implements a Runge-Kutta technique with a variable 
time step for efficient computation. However, this was preceded by the assessment of the 
values of R, and R, at œ= 0.00001, o =6 ,c=10 , d =0.00014, u =0.486, ¢ =0.96 and 


q=2.We found that0.2626= R, < R, =1.3129, implying that therapy reduces the HCV 


infection in vivo. 


According to Fig. 2, the simulation for the basic reproduction number R, shows that the 
HCV transmission is very high without therapy, regardless of variations of spontaneous 
recovery of infected hepatocytesg, for0<q<2.In Fig.3, the numerical solution of R, 
shows four different cases: (1) when q =Oand ¢=0, R,has the highest value, indicating 
the increase in the chance of disease evolution to chronic state, implying that patients who 
do not spontaneously clear the virus and are not subjected to treatment will get chronic 
HCV infection in the long run; (2) when q =Oand ¢+0,R, is still high although decreases 
with increasing value of € ,implying that patients who do not spontaneously clear the virus 
will only clear it through drug therapy. But, the drug will totally eradicate the disease with 
increase of blocking efficacy of the virus production; (3) when q #0ande =0,R,is 
relatively low, implying that patients who spontaneously clear the virus and are not 
subjected to therapy will cure from the disease at different times depending on the rate of 
clearance qand (4) wheng #0and ¢+0,R,attains the lowest value and additionally 
decreases with increasing values of q and €, implying that patients who spontaneously 
clear the virus and are subjected to therapy will also cure from the disease. The numerical 
solution based on category 4 shows a very low value of R,, implying that the disease is 


eradicated very rapidly. However, this is practically not recommended as therapy of this 
group of patients is unnecessary due to spontaneous clearance of the virus they undergo. 
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Figure 2: Relative values of reproduction numbers R and R, with respect to the rate of 


spontaneous recovery of infected hepatocytes (q). 
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Figure 3: Value of R, with respect to drug blocking efficacy of the virus production rate (e) , 


with variation in the rate of spontaneous cure of infected hepatocytes (q =0,1,2 ) 












































= o8 

3 — q = 0 
E — g = 
2 0.6 = 
£ 

o 

8 04 

o 

S 

2i Ra 

E 

E 5 1 2 3 4 5 6 


Rate of virus production(c) 


Figure 4: Value of R,with respect to the rate of virus production rate (o) , with variation in 


the rate of spontaneous cure of infected hepatocytes (q =0,1,2) at £ = 0.96 
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Figure 5: Value of R, with respect to rate of virus production rate (o) , with variation in the 


rate of spontaneous cure of infected hepatocytes (q =0,1,2) at £ = 0.81 


IMMUNE SYSTEM RESPONSE AND THERAPY 19 


e 









































Effective reproductive number(R) 





o 1 2 3 4 5 6 


Rate of virus production(c) 


Figure 6: Value of R, with respect to the rate of virus production rate (co) , with variation in 


the rate of spontaneous cure of infected hepatocytes (q = 0,1,2 ) at £ =0 


Fig.4, Fig.5 and Fig.6 show that the value of R, increases with the rate of virus production 


oate =0.96,¢=0.8land ¢ =Orespectively, and it attains the lowest and greatest values at 
q=2and q =0 respectively. It also reduces with increasing value of ¢.This implies that 


patients with different rates of viral production and who do not spontaneously clear the 
virus can remove it by means of only therapy; this is optimal if the dug blocking efficacy is 
approximately 100% (e.g £ = 96% ). However patients with higher rates of viral production 
(e.g.o0 =2) cure later than those with lower rates (e.g.o =1).Conversely, patients who 
spontaneously clear the virus will do it faster with therapy. It will vanish in a short time if 
the drug is highly effective. Practically, therapy is necessary and so should be directed to 
only patients who do not spontaneously clear the virus. 


6. Conclusion 

In this paper, we have formulated a deterministic mathematical model on the dynamics of 
HCV and immune system with the effect of therapy in the acute phase of infection. We have 
shown that the numerical solution of the basic and effective reproductive number obtained 
validate that therapy has significant impact on the disease transmission, which implies that 
therapy reduces the number of secondary infections produced by an infectious hepatocyte 
during its infectious period, i.e. fewer hepatocytes are infected. From the sensitivity 
analysis, the parameter for the drug blocking efficacy is the most negatively sensitive 
parameter while the parameters for the rates of infections, virus production and 
recruitment of susceptible hepatocytes and mortality rate due to infection are the most 
positively sensitive ones. Thus, if the blocking drug efficacy of the virus production 
increases there is reduction of disease transmission, and eventually it dies out. 


Literatures reveal that therapy of acute HCV infection is associated with challenges such as 
unclear criteria of identification of patients who require therapy (to exclude those who are 
able to spontaneously clear the virus), initiation of therapy and therapy regimens. 
Available drugs for therapy of chronic infection have severe side effects due to combination 
therapies. Nevertheless, therapy of people with acute HCV infection is necessary to prevent 
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evolution of the disease to chronicity, which is a fatal period of infection. Current models 
describing the dynamics of HCV with therapy do not consider drug types, variation of 
blocking efficacy with dosage and the virus genotype concurrently. For example, 
individuals with genotype 1 and 4 HCV infections do not respond to most antivirals while 
those infected with genotype 2 and 3 HCV respond quite well to these antivirals, and 
therefore the mathematical models are only descriptive of the dynamics of HCV with 
therapy for the second category of patients. Hence, more mathematical models are needed 
to address all these anomalies concurrently. However, our model assumes that the drug 
blocking efficacy of virus production increases with drug dosage administered to a patient 
who is a responder to the drug and can either clear the virus spontaneously or not. 


Hence we recommend that more mathematical models should be formulated to describe 
viral dynamic profiles observed when new antiviral drugs are employed, test possible 
mechanisms for new biological or clinical phenomena and optimize therapy such as 
determining optimal drug dosages and therapy duration. The models should be used to 
give additional insight into the HCV pathogenesis and assist in the development of more 
effective therapeutic strategies. The virus genotype should be considered when modeling 
the HCV dynamics and immune system with therapy, as the current mathematical models 
assume that all patients are good responders to the available antivirals. Also, they should 
consider the dynamics of HCV in a patient with a graft organ such as kidney, liver or heart. 
Modeling the dynamics of HCV and immune system with drug therapy can be done to 
include co-infections such as HCV-HIV co-infection. 
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